Universalities in ultracold reactions of alkali polar molecules 
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We consider ultracold collisions of ground-state, heteronuclear alkali dimers that are susceptible 
to four-center chemical reactions 2 AB — > A2 + B2 even at sub-microKelvin temperature. These 
reactions depend strongly on species, temperature, electric field, and confinement in an optical 
lattice. We calculate ab initio van der Walls coefficients for these interactions, and use a quantum 
formalism to study the scattering properties of such molecules under an external electric field and 
optical lattice. We also apply a quantum threshold model to explore the dependence of reaction 
rates on the various parameters. We find that, among the heteronuclear alkali fermionic species, 
LiNa is the least reactive, whereas LiCs is the most reactive. For the bosonic species, LiK is the 
most reactive in zero field, but all species considered LiNa, LiK, LiRb, LiCs, and KRb share a 
universal reaction rate once a sufficiently high electric field is applied. 
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I. INTRODUCTION 

The study of ultracold polar molecules has now be- 
come a vast and exciting area of interest since the forma- 
tion of bi-alkali heteronuclear polar molecules [U-IE] . The 
molecules can be controlled at the ground electronic, vi- 
brational, rotational and hyperfine @ quantum-level. 
The external motion of the polar molecules can also be 
modified by an electric field [7| and by an optical lattice 
confinement Q. 

Polar molecules offer remarkable characteristics. First, 
they have strong electric dipole moments [l(|. The 
interactions between polar molecules can then be domi- 
nated by electric dipole-dipole terms. The electric molec- 
ular interactions are strong, long-range, anisotropic and 
can be tuned by electric fields. Secondly, the polar 
molecules can be either bosons or fermions. If the polar 
molecules are addressed in a single quantum state, they 
become indistinguishable and quantum statistics plays 
a strong role. An ultracold gas of bosonic molecules 
can lead to Bosc-Einstein condensation and an ultracold 
sample of fermionic molecules can lead to a Degenerate 
Fermi gas . Thirdly, two polar molecules can be reactive 
or not [lll4l^ | . It was found in Ref . [ll| that among the 
bi-alkali heteronuclear molecules in their absolute funda- 
mental ground state, that the Lithium species LiNa, LiK, 
LiRb, LiCs in addition with the KRb molecule (category 
1) gave rise to two-body exoergic chemical reactive pro- 
cesses while the remaining species NaK, NaRb, NaCs, 
KCs, RbCs (category 2) resulted in two-body endoergic 
processes. Reactivity is an adva ntag e to investigate the 
ultracold chemistry of molecules [15[ ■ It also provides a 
clear signature (in term of molecular loss) of two-body 



interactions in a gas and depends strongly on the ap- 
plied electric field [161 ]. The non-reactive molecules have 
the advantage of being chemically stable in their absolute 
ground state and can help to reach long-lived samples of 
polar molecules. However, if dense samples of molecules 
are formed in Bose-Einstein condensates for example, 
three-body collision can become a source of loss, and it 
is important to investigate the collisional properties of 
such processes p7HT9l ] . Finally, molecules offer a rich in- 
ternal quantum structure and can be manipulated with 
electromagnetic waves in order to address their quan- 
tum state. Exciting perspectives have been proposed for 
these polar molecules. This involves condensed matter 
and many-body physics, quantum magnetism, precision 
measurements, controlled chemistry and quantum infor- 
mation (13-0. 
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For all these reasons, many experimental groups are 
currently interested in creating polar molecules. The 
fermionic polar molecules 40 K 87 Rb received a particu- 
lar experimental @, M EE E3 and theoretical [U, [H- 
[IH consideration recently. However, much less is known 
about the interactions and the dynamical properties of 
the other polar bi-alkali molecules, for which experimen- 
tal attention is also devoted [1 @, H HiHIi]. This is 
what we address in this article. In Section II, we com- 
pute the isotropic long-range van der Waals coefficients 
between polar molecules. We focus our study to the ex- 
oergic molecules (category 1). In Section III, we use 
these parameters to perform quantum scattering calcu- 
lations assuming full loss when the polar molecules are 
close to each other. We consider the case of collisions 
in free and confined space, in electric fields. We use a 
Quantum Threshold (QT) model to explain how the col- 
lisional properties scale with the different species. We 
arrive at analytical expressions of high-loss collision rates 
of bosonic or fermionic molecules, which can also be ap- 
plied to the inelastic and reactive case of molecules of 
category 2, as well as atom-atom or atom- molecule colli- 
sions, provided the van der Waals coefficients arc known. 
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We conclude in Section IV. 
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II. ISOTROPIC LONG-RANGE INTERACTION 
OF REACTIVE POLAR MOLECULES 

The isotropic dispersion coefficient Cq between two 
identical diatomic alkali-metal molecules in the v=0 and 
J=0 rovibrational ground state of the X 1 S + potential 
has three contributions 
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where the first term of the integrand is the square of 
the isotropic dynamic polarizability a gr (iuj) at imaginary 
frequency ito due to rovibrational transitions within the 
ground state potential. The second term in the integrand 
is the square of the isotropic polarizability a exc (ioj) due 
to transitions to the rovibrational levels of electronically 
excited potentials, while the last term indicates an in- 
terference between the first two contributions. In these 
and subsequent expressions both the dispersion coeffi- 
cient and the polarizability are in atomic units. A thor- 
ough discussion of dispersion forces between molecules 
can be found in Ref. |47| . 

We find that a gx (iu)) = a 0g /(l + (uj/r] g ) 2 ) [H EH to 
good approximation with a^ g = d 2 /('SB) and j] g = 2B, 
where d p and B are the electric permanent dipole mo- 
ment and rotational constant at the equilibrium sepa- 
ration R e between the atoms in the molecule, respec- 
tively. The contribution from transitions between vibra- 
tional levels within the ground state potential is negligi- 
bly small. Consequently, = dp/(6B) in agreement 
with the findings of Ref. J49| ■ 

The isotropic dynamic polarizability a 0KC (iuj) contains 
contributions from transitions to the rovibrational ex- 
cited 1 E + and 1 II potentials, which correspond to the 
parallel and perpendicular component of the polarizabil- 
ity, respectively. Based on the Franck- Condon prin- 
ciple we can evaluate the polarizability at each inter- 
atomic separation rather than perform an average over 
ro-vibrational levels [Efj. For v=0 and J=0 the sep- 
aration is R = R e . We then parametrize a exc (iuj) ~ 
J2j a j with 
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Each term corresponds to an excited potential. In prac- 
tice, we have found it more convenient to evaluate the 
polarizability at R = R e as function of real frequen- 
cies and find the parameters aoj and rjj from a fit. The 
static polarization due to the excited state potentials is 
Q!exc(0) = J2j a 0j- Using Eqs. (fT]) and ((2]) we obtain the 
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The dynamic polarizability a C xc(w) at real frequency 
is calculated using a coupled cluster method with sin- 
gle and double excitations (ccsd) [5l[. The calculation 
of the static polarizability and permanent dipole mo- 
ment is performed at much higher level using coupled 
cluster method with the single, double and triple excita- 
tions (ccsdt). Twelve electrons, including ls 2 2s 1 of the 
Li atom and (n— l)s 2 (n— l)p 6 ns 1 of the Na, K, Rb, and 
Cs atoms, were explicitly used in both ccsd and ccsdt 
calculations. The dipole moment for each molecule was 
averaged on the zero vibrational level. We employed the 
cc-pCVQZ basis sets for Li and Na from Refs. Il52l . l53j , 
the all-electron basis for the K atom from Ref. |54( , and 
the ECP28MDF and ECP46MDF basis sets with the rel- 
ativistic effective core potentials from Ref. (5|| for the 
Rb and Cs atoms. A comparison of our data on the 
dipole moment and static polarizability with results of 
Refs. shows a good agreement within a few %. 

Table |T] lists our Cq coefficients for four pairs of iden- 
tical alkali-metal molecules in the v = 0, J = rovibra- 
tional level of the X 1 S + potential. For completeness, we 
tabulate the contribution to the isotropic component of 
the static polarizability from electronically excited poten- 
tials, the rotational constant, and the permanent dipole 
moment for each of the four molecules. 

Tabic [I] shows that the value of the Cq coefficient as 
well as the three contributions to it increase when we 
move down along the first column of the periodic table 
for the second atom in our four diatomic molecules. Most 
of the increase can be traced back to increasing perma- 
nent and transition dipole moments. For example, the 
ground state contribution Cg SI ' increase by four order of 
magnitude as the permanent dipole moment increases by 
a factor of ten. For the excited state contribution Cg exe -* 
the increase is less dramatic as the transition dipole mo- 
ments increase only weakly. Only for the LiNa molecule 
does the excited state contribution dominate the Cq co- 
efficient. 



III. DYNAMICS IN THREE-DIMENSIONAL 
SPACE 

A. Quantum numerical calculation 

We use the isotropic van der waals Cq coefficients cal- 
culated in the previous section to compute the chemi- 
cal rate coefficients of the reactive polar molecules. Wc 
use the same formalism used in Ref. Q for the chem- 
ical reaction KRb + KRb — » I<2 + Rb2. Wc employ 



TABLE I: Van der Waals Cq coefficients in atomic units for 
the interaction between the two molecules in the v — 0, J = 
rovibrational levels of the X 1 S + potential and other molecu- 
lar characteristics used to calculate Cq \ a e xc(0) is the isotropic 
static polarizability due to transitions to electronically excited 
potentials; B and d p are the rotational constant and electric 
permanent dipole moment, respectively. These three prop- 
erties are evaluated at the equilibrium separation R e . The 
value of B is from [13] . The next three columns are the ex- 
cited state, interference, and ground state contributions to 
the total C*6, shown in the last column. 
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FIG. 1: (Color online) Solid lines: Loss rate coefficient /3l=i 
divided by T in free 3D space, of the reaction AB + AB — > 
A-2 + B2 for different reactive fermionic polar molecules AB 
= LiNa, KRb, LiK, LiRb, LiCs, as a function of the elec- 
tric dipole moment. The fermions are considered in a same 
indistinguishable quantum state. Dashed lines: /3l=i,m l =o 
and Pl=i,\m l \=i components of L — 1, shown here for AB = 
LiNa. 



where L is the quantum number associated with the or- 
bital angular momentum of the collision, and Ml, the 
quantum number associated with its projection onto a 
quantization axis (see Ref. for details). Solving the 
eigenstates of the Hamiltonian leads to the set of close- 
coupling equations 



a time-independent quantum formalism, including only 
one molecule-molecule channel corresponding to the ini- 
tial state of the molecules, but including several partial 
waves. For two particles of mass to 1,771,2, the Hamilto- 
nian of the system is given by 



H = T + V ahs + Vy 



dW 



Kid- 
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Using spherical coordinates (r, 9, </>), the kinetic energy is 
T = — fi 2 V^/(2/i) , fi = toiTO2/(toi +TO2) is the reduced 

mass of the colliding system, V a bs = iAe~^ r ~ rmi "^ rc is 
an absorbing potential to account for the loss of particles 
due to chemical reactions or inelastic collisions in the in- 
cident channel, where A is the strength of the absorbing 
potential, r min is the position where the potential starts, 
and r c is the position where the potential vanishes expo- 
nentially. Vvdw = —Cq/R 6 is an isotropic van der Waals 
interaction, and Vdd = [^1^2 (1 — 3 cos 2 8)}/(4Tre R 3 ) is 
the dipole-dipole interaction between the two particles if 
an electric field is applied. Here di,d,2 are the induced 
electric dipole moments in the laboratory frame and their 
maximum value is given by their permanent dipole mo- 
ment cJ Pl i,d Pj 2 m the molecular frame. We expand the 
total wavefunction onto a basis set of spherical harmonics 
(or partial waves) 

^ (R, 0^) = ^J2 Y l' L ( R > °) F l> L ( 6 ) 
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E represents the total energy which is, in this study, the 
collision energy E c , as we use only one molecule-molecule 
incident channel. We use the same notation as in Ref. [l6| 
C 3 (L, L>\ M L ) = a0, d 1 d 2 /47re with 
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The effective potential in Eq. ([7]) is given by 

h 2 L(L + l) C 6 C 3 (L,L;M L ) 



V cff 



2nR 2 



R 6 



R 3 



(8) 



(9) 



for a given L,Ml- The absorbing potential is chosen 
in Eq. (|7J) in such a way that the elastic probability 
vanishes (or the loss probability is unity) when the two 
molecules come close together. The case for which the 
loss probab ility is smaller than unity has been discussed 
in Ref. [II il 11 . 
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the LiCs system. Now the rates seem to scale with the 
reduced mass of the system only. For a given dipole, 
the electric dipole interaction is the same between the 
species, only the centrifugal terms differ. Higher mass 
means smaller barrier so higher loss rate. 

For the bosonic case, in the van der Waals regime, 
KRb are the least reactive molecules, followed by LiNa, 
LiRb, LiCs and finally LiK. Bosonic particles collide in 
a s-wave at ultralow energy where no incident barrier 
is present. Instead, one must invoke the probability for 
quantum transmission. In the electric field regime, all 
different systems have the same rate coefficients. This 
will be explained in the next section. 



FIG. 2: (Color online) Same as Figure[T]for different reactive 
bosonic polar molecules. The rate coefficient /3l=o is plotted 
as a function of the electric dipole moment for T — > (Wigner 
regime). The bosons are considered in a same indistinguish- 
able quantum state. 



We report in Fig. Q] and [2] the loss rate coefficient as 
a function of the induced electric dipole moment d, for 
two indistinguishable fcrmionic molecules (Fig. [T]) and 
for two indistinguishable bosonic molecules (Fig. [2]), for 
LiNa-LiNa, KRb-KRb, LiK-LiK, LiRb-LiRb and LiCs- 
LiCs collisions. To converge the results, we use five par- 
tial waves, L = 1,3,5,7,9 for the fermions and L = 
0, 2, 4, 6, 8 for the bosons. We used the values of Cg and 
dp reported in Tab. HI and the value of Cq = 16133 a.u. 
of Ref. [U and d p = 0.566 D of Ref. j| for KRb. We 
provide a list of the fermionic and bosonic isotopes of 
each species in Appendix A. These results have been ob- 
tained in the regime of ultracold temperature. In this 
regime, the fermionic rate scales linearly with the tem- 
perature (hence we have plotted the rate divided by the 
temperature) while the bosonic rate is independent of the 
temperature according to the Bethe- Wigner laws [56|, [57| ■ 
For both cases, the rate scales as a constant in the van 
der Waals regime where d —¥ 0, and an increasing term 
in the electric field regime where d — > d p . We note that 
for large dipole moments, the corresponding dipole length 
«dd = fJ,d 2 /h 2 may exceed the distance between molecules 
given by the inverse third of the molecular gas density 
<ham = nr 1 / 3 . In such situation, there are no more col- 
lisions between molecules. Instead, a dense liquid/solid 
phase is entered where many-body physics becomes im- 
portant. 

For the fermionic case, in the van der Waals regime, it 
is seen that the LiNa system is the least reactive, followed 
by KRb, LiK, LiRb and finally LiCs. Qualitativelly, light 
masses and small values of Cg increase the incident p- 
wave barrier (this is the case for LiNa) and hence decrease 
the chance to get high chemical reactivity, while heavy 
masses and large values of C§ decrease the barrier (this 
is the case for LiCs) and increase the reactivity. In the 
electric field regime, the same general trend is observed, 
except now the rate of the KRb system is as high as 



B. Quantum Threshold model 

To understand the physical trends seen in the numeri- 
cal results, we employ an analytical Quantum Threshold 
model (QT model) [161 ] which provides a universal ex- 
pression of an ultracold collision (chemical reaction or 
inelastic collision) with short range unit loss probability. 
The QT model is a clear and simple model to describe the 
dependence of an ultracold chemical reaction on the re- 
duced mass and the isotropic van der waals Cq coefficient 
of the molecule-molecule complex, and on the induced 
dipole moment via the presence of an applied electric 
field. The QT model assumes that the loss probability 
scales as 
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Pl.IMlI 1 E ^ 
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(10) 



where E* is a characteristic energy corresponding to the 
long range interaction of the molecules in a partial wave 
L,Ml. p L |M ! is a dimensionless quantity of order of 
unity, and is estimated by fitting the expression with the 
numerical results. The thcrmalizcd rate coefficient is ex- 
pressed by 



L,M L = P L 




(11) 



where the brackets denote a Maxwell-Boltzmann distri- 
bution over the collision energy to the power L. A = 2 
if the particles are in indistinguishable states and A = 1 
if they are in distinguishable states [58| . 



1. QT model for p-wave collisions 

For p-wave collision (L = 1), we chose the characteris- 
tic energy E* equal to the height of the incident barrier, 
, j 01 the effective potential V e g, composed of the 
strongest attractive potential —C n /R n and the strongest 
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repulsive potential C m /R 

( 



mC„ 



a 



y rn 



(nC n \ 
\mC m ) 



(12) 



The position of the barrier is given by 

nC n 



m C„ 



(13) 



The combinations of n and m are given in Tab. |TT] with 
the corresponding height of the barriers. For the van 
der Waals regime and for either \Ml\ = 0, 1, the height 
of the barrier is made by the the attractive van der 
Waals interaction — Crj/i? 6 and the repulsive centrifugal 
term C 2 /R 2 = h 2 L(L + l)/(2/xR 2 ) with L = 1, giv- 

6 2 

ing rise to a characteristic energy 1 y For the elec- 
tric field regime and for |Ml| = 0, the height of the 
barrier is made by the attractive dipole-dipole interac- 
tion -C 3 (l,l;0)/i? 3 = -[(4/5)rf 2 /(47T£ )]/i? 3 and the 
repulsive centrifugal term C2/R 2 = h 2 L(L + l)/(2^i? 2 ) 
with L = 1, giving rise to a characteristic energy -E^'q. 
Finally, for the electric field regime and for \Mj,\ = 
1, the height of the barrier is made by an attrac- 
tive -C A /R i = -[(72^/(875ft 2 ))d 4 /(47re ) 2 ]/i? 4 and 
the repulsive dipole-dipole interaction —(73(1, 1; 1)/-R 3 = 
+ [(2/5) d 2 / (Atteq)]/ R 3 , giving rise to a characteristic en- 



The —Ci/R attractive interaction comes 



er §y £-1,1- 

from the coupling between the L = 1 and L = 3 of the 
\M h \ = 1 component. This is demonstrated in Appendix 
B. 



TABLE II: Characteristic energies E r for L = 1, \M L \ = 0, 1 
in the van der Waals (vdW) and electric (elec.) regime. 
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\Ml\ 
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C m /R m 
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L = 1.\M L \ 
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25 ffi ( d 2 \~ 2 
108 M a V ilTe ° ) 
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4-n-eo R 3 
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Replacing these three values of E* into Eq. (fTT|l for 
L = 1 and assuming (_E C ) = 3/cbT/2, where fcs is the 
Boltzmann constant and T the temperature, we arrive 
at the following expressions for the \Ml\ = 0, 1 rate as 
d —> in the van der Waals regime 
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FIG. 3: (Color online) Top panel: Van der Waals regime for 
Ml = and \Ml\ = 1. The quantity = divided by 
T is plotted as a function of (/x Ce) 3 ^ 4 . Middle panel: Electric 
field regime for Ml = 0. The quantity f3f% divided by d 6 T 
is plotted as a function of /i 3 . Bottom panel: Electric field 
regime for \M L I = 1. The quantity /3f'f divided d 6 T is plotted 
as a function of fx . 



(15) The I Mi I = rate as d —> d p in the electric field regime 
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o elec 



3tt /6 
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k B TxA (16) 



with 



Pi'o = 0.54 ±0.04. 



(17) 



Finally, the | Ml | = 1 rate as d — > d p in the electric field 
regime is given by 



o elec 

J £=l,|Mi,|=l 



p 4 ' 3 — (20000(24/875) 
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with 



p\'\ = 0.16 ±0.02. 
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The coefficients p"^ r ^ Mh 



energies E n,m 
results in Eq. 



associated with the characteristic 
are found by confronting the analytical 
HI), Eq. (HU), and Eq. (UHl) with our nu- 



merical calculations of Fig. Q] The quantity Pl=i\m l \=q i 
divided by T obtained from the numerical results, is plot- 
ted as a function of the quantity (/i C§) 3 / 4 for the van der 
Waals regime in the top panel of Fig. 03 The quantities 

^=1\ Ml \=o and P£=t\M L \=i dividcd b y d& and T are 

plotted as a function of the quantity [i 3 for the electric 
field regime for the \Ml = 0| and \Ml = 1| component 
in the middle and bottom panels of Fig. [3] respectively, 
for the different fermionic reactive systems. We find that 
the numerical results fit a line, confirming the validity 
of the QT model analysis (the fitting uncertainty of the 
lines provides an uncertainty to the P L ' \m l \ parameters). 
The fitting parameters are the slope of these lines and 
are reported in Eq. ([15]), Eq. ((I7|) and Eq. (fT9]) . 

We see that both components \Ml\ =0,1 analytical 
rates (Eq.([T4])) at ultracold temperature are the same in 
the van der Waals regime and are dictated by a d 6 depen- 
dence in the electric regime (Eq. (fl6|) and Eq. (fT8|) ) with 
different magnitudes. These expressions provide a clear 
explanation of the trends observed numerically. The loss 
rate behaves as (/xCg) 3 / 4 in the van der Waals regime. 
In the electric field regime, the loss rate scales as /i 3 , 
increasing only with the mass. In both regimes, these 
expressions explain why fermionic LiNa is the least reac- 
tive alkali polar species and fermionic LiCs is the most 
reactive one. 

0.53 ± 0.07 is 



6.2 



(0,1) 



We note that the results of p l 
in very good agreement with the analytical expression of 
2 19 / 4 tt/(3 17 / 4 [r(3/4)] 2 ) = 0.528 found using a Quantum 
Defect Theory (QDT) The values p 3 { 2 = 0.54 ± 

0.04 and p 4 ' 3 = 0.16 ± 0.02 for the 1/i? 3 interaction in 



the electric field regime have not to our knowledge been 
determined analytically in a QDT framework. 

We also note that these constants barely change be- 
tween the regime dominated by the van der Waals in- 
teraction and the regime dominated by an electric field 
interaction for the \Ml\ = component. The ratio of the 
\Ml\ = 1 over the |Ml| = component in the electric 
field regime is 0.003. As a consequence, the \Ml\ = 1 
component is negligible in the electric regime, as seen in 
Fig. [Q for the LiNa system, and one can provide an esti- 
mation of the total p-wave rate coefficient for the reactive 
systems by 

Pl=i = Pl=i,\m l \=v> + ^ Pl=i,\m l \=i 

vdW , a el 



o a vdW 1 o elec 

o Pi=l,|Af I .|=0,l "T" Pl=1,|M l |=0 

7T 



0.53 
+ 0.54 x 



3 17 M 3 a 3xl/4 



h 10 

/ 2 9 / 2 3 n / 2 fi 3 

\ 5 3 h 7 



d 6 



(4^ ) 3 



k B T x £20) 



2. QT model for s-wave collisions 

For s-wave collisions (L = 0, Ml = 0), there is no inci- 
dent barrier because the repulsive centrifugal term van- 
ishes. It is possible however to estimate a characteristic 
length and energy [59[ given respectively by 



2^C„ 

— w 



e': 



'(tat, 



(21) 



In the van der Waals regime, the characteristic energy 
is E% = h 3 /y / 2 3 (i 3 C 6 - In the electric field regime, the 
electric dipole-dipolc interaction vanishes for L = 0. But 
as there is a coupling between the L = and the L — 2 
component in Eq. ([5]), it is found after diagonalisation, 
that the electric dipole interaction behaves as a —C4/R 
with C 4 = 4^d 4 /[15?i 2 (47T£ ) 2 ] (sec Appendix C). In 
return, this corresponds to a characteristic energy Eq = 
15 h 6 (47T£ ) 2 /[16 [i 3 d 4 }. This is summarized in Tab [EH] 



TABLE III: Characteristic energies E, for L = 0, \M L \ = in 
the van der Waals (vdW) and electric (elec.) regime. 



regime 


Ml 


- C n /R n 


l=o,\Mt. 1 


vdW 
elec. 






c 6 
" W 

4 /id 4 


V 23 c 6 

15h 6 ( d 2 \~ 2 


15 h 2 (4ire ) 2 -R 4 


16 M 3 \4,Ttea J 



Replacing these two values of E* into Eq.lfTTjl for 
L = 0, we arrive at the following expression for the 
L = 0, \Ml\ = rate as d — > in the van der Waals 
regime 



8 vdW - n 6 7T ( 



1/4 



x A 



(22) 
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FIG. 4: (Color online) Top panel: Van der Waals regime. The 
quantity /3o,o is plotted as a function of (Ce/n 3 ) 1 ^ 4 - Bottom 
panel: Electric field regime. The quantity /3q divided by d 2 
is plotted for the five different colliding species. 



with 



pl o = 1.92 ±0.01. 



(23) 



The L = 0, | Mi | = rate as d — > d p in the electric field 
regime is 



^clcc 4 w VW30 d 2 



h 47T£ 







with 



Po.o 



3.74 ±0.05. 



(25) 



Compared to L — 1, the rates at ultracold temperature 
for L = behave now as (Cq/ p 3 ) 1 ^ 4 in the van der Waals 
regime, making bosonic KRb molecules the least reac- 
tive ones and bosonic LiK molecules the most reactive 
ones, due to the interplay between the Cq coefficients 
and the cube of the mass. In the electric field regime, 
the rates behave as d 2 and are independent of the mass, 
so that for the same induced dipole, all bosonic polar 
molecules react with the same rate coefficient. The coef- 
ficients j»2 I Mi, I associated with the characteristic energies 

E 1,\m l \' are found b y P lottin g the quantity ^l=q!\m l \=o 



obtained from the numerical results of Fig. [2] as a func- 
tion of the quantity (Cq/ '/i 3 ) 1 / 4 for the van der Waals 
regime in the top panel of Fig. [H and the quantity 
P£=o \m l \=o divided by d 2 for the electric regime in the 
bottom panel of Fig. [4] for the different bosonic reactive 
systems. As for the fermionic case, the numerical results 
form a line for the first plot and are constant for the sec- 
ond plot, validating the QT model analysis. Again, we 
note that the results of pi = 1.92 ± 0.01 in Eq. ((23j) 
is in very good agreement with the analytical expression 
of 8 7r/[r(l/4)1 2 = 1.912 found using a Quantum Defect 
Theory [29| or a Quantum Langevin Theory (QL) [3o| . 
The value of Pq = 3.74 ± 0.05 agrees within 7% with 
the analytical expression of 4 from a Quantum Langevin 
Theory [60( using the —C4/R 4 interaction in the electric 
field regime. One can formulate a good approximation 
for the s-wave loss rate coefficients by 



o a vdW j_ a cloc 

PL=Q — Pl=o,\m l \=o ~t Pl=o,\m l \=o 



clcc 



TT< 1.92 X 



3.74 x 



2h 2 C 6 

VWso d 2 



1/4 



47T£r 



x A. 



(26) 



The formulas from Eq. {14]) to Eq. ([19]) and from 
Eq. (f22j) to Eq. (|25|) can be used to determine the 
inelastic and reactive collisional properties of other 
atom-atom, atom-molecule or molecule-molecule colli- 
sions, provided that full loss occurs when they encounter 
one another. This case can occur for molecules of 
category 2 (NaK, NaRb, NaCs, KCs, RbCs) if the 
molecules are not in their absolute ground state, for 
example in a higher vibrational state, where inelastic 
molecule-molecule collision can occur or when the 
reactants have higher energy than the products so 
that an exoergic reaction can take place. What is 
left unknown is the Cq coefficients (except for RbCs), 
for each of these initial ro-vibrational states of these 
molecules and has to be calculated individually. For the 
RbCs molecule the Cq coefficients have been calculated 



as a function of vibrational quantum number in Ref. [31 



We provide in Appendix D the corresponding QT ex- 
pressions for the imaginary part of the scattering lengths 
for s-wave collisions and scattering volumes for p-wave 
collisions. 



IV. DYNAMICS IN TWO-DIMENSIONAL 
SPACE 

For the confined 2D scattering we use the same formal- 
ism developed in Ref. [HI, [36| . The confinement is given 
by an optical lattice in the z direction, which we approx- 
imate by a harmonic oscillator potential Vho = fiuj 2 z /2 
of frequency v and angular frequency uj = 2 ir v. One 
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FIG. 5: (Color online) Loss rate coefficients for different reac- 
tive polar molecules in confined 2D space, for indistinguish- 
able fermions (top panel) and indistinguishable bosons (bot- 
tom panel). The frequency of the ID trap is v = 20 kHz 
and the temperature is T = 500 nK. The dashed lines rep- 
resent a model based on the rescaled 3D rate coefficients for 
add "C aho, and the dotted lines represent a model based on 
a functional form (Ref. [38| for example) for add 3> flho- 



can also define a harmonic oscillator confinement length 
aho = yh/JJiuj). We consider the dynamics of two 
molecules in the ground state of this harmonic oscillator. 
In confined space, Ml remains a good quantum number. 
Additional selection rules apply and for indistinguishable 
bosons, \Mj\ = 0, while for indistinguishable fermions, 
| Ml | = 1 [111, [lH, for molecules in the ground state of 
the harmonic confinement. We present in Fig. [5] the loss 
rate coefficient for a confinement of v = 20 kHz as a 
function of the dipole moment for a given temperature 
T = 500 nK, for the fermionic species (top panel) and 
the bosonic species (bottom panel). We use fourty par- 
tial waves, L = 1 — 79 for the fermions and L = — 78 
for the bosons, to converge the results. At small elec- 
tric dipoles, when add <C <Xho, the collisions are quasi-2D 
(q2D) and the loss rate coefficients display a similar be- 
havior than their 3D counterpart for |Mt| = 1 for the 
indistinguishable fermions and for \Ml\ = for the in- 
distinguishable bosons. At large electric dipoles and for 



LiK, LiRb, LiCs, when add 3> aho, the collisions are fully 
2D and the loss rate coefficients show a suppression as 
discussed in Ref. [U [HQ . 

For the quasi-2D regime add *C aho, we compare in 
dashed lines in Fig. [5] a two dimensional loss rate coeffi- 
cient rescaled from the numerical calculation in three di- 
mensions [HI, Hi], HH from the previous section for L = 1 



3 B 3D 

pq2D _ £ A 



j3D 



2 A/n^ho 



(27) 



where the factor 3/2 accounts for the difference of the 
mean energies in 3D and 2D for a given temperature T 
(in 3D, (E c ) = 3k B T/2 while in 2D, (E c ) = k B T). For 
L = 0, we get 



P 



q2D 



/7T a h( 



7T H 



(28) 



We found for the fermions a good agreement between 
the numerical 2D rates (in solid lines) and the rescaled 
from 3D rates (dashed lines). For the bosons, a good 
agreement is found for the LiNa system, but not for the 
other systems like KRb for example, even if the order 
of magnitude is right. For bosons, threshold laws dis- 
play a logarithmic dependence and are not accounted in 
Eq. (|28p . The QT formulas which describe the numeri- 
cal 3D rates can also be rescaled in the same manner so 
that a good approximation to the loss rate coefficient for 
fermions in the quasi-2D regime add <C aho is given by 



P L =1 - 2X n\^j;P\M 




0.16 x 



d 6 



3tt / 20000 72 



8 V 27 875 / 



k B T x A. 



3/2 3 

h 7 



For bosons, the rescaled QT formula is 



^L=0 — 



(29) 



n riuj \ I 2 h Cfi 

' 1.92 x ' 



1/4 



3.74 x 



x A 



(30) 



but is only a good approximation for LiNa. 

For the 2D regime add &ho, we compare in Fig. [5] a 
functional form provided in Rcfs. [3(3, [H, H3, f38l |63H65| . 
We found that the forms 



f3 2D = 2 x 13 - 



[i \huj/2 



-2(a dd /a ho ) 2 / 5 x A ( 31 ) 
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for indistinguishable fcrmions, and 



Acknowledgments 



j2D 



13- 



H \hu)/2 



-2(a dd /a ho ) 2 / 5 



X A 



(32) 



for indistinguishable bosons fit well the numerical data. 
These formulas are reported in dotted lines in Fig. [5] 
We find a coefficient of f3 in front of the exponential 
and a coefficient of 2 inside the exponential, by fitting 
our numerical results. These values are different from 
the values found in Ref. [38|. This is attributed to the 
different regimes of collision energies and confinements 
involved in the fitting. It has been shown in Ref. (38| 
that the fitting parameters of the functional form may 
differ for different values of the collision energies. 



V. CONCLUSION 

By computing the Cq coefficients for different pairs 
of alkali polar molecules of LiNa, LiK, LiRb, and LiCs, 
and using an available one for KRb, we estimated the 
quenching rate coefficient assuming full loss when they 
encounter one another, for the fcrmionic species and for 
the bosonic species, both for the van der Waals regime 
and the electric field regime. We found that, at ultra- 
cold temperature, fermionic LiNa is the least reactive 
system while LiCs is the most in the van der Waals 
regime and electric field regime, due mainly to the in- 
crease of the Ce coefficient for the former regime and 
due to the increase of the mass for the later. Bosonic 
KRb molecules are found to be the least reactive ones 
while LiK the most in the van der Waals regime. All the 
bosonic molecules were found to have the same universal 
reactive rate in the electric field regime. These behaviors 
were all explained using a Quantum Threshold model. 
From our numerical results, we found analytical expres- 
sions for the reactive rate coefficients for fermionic and 
bosonic molecules, in the van der Waals and electric field 
regime. These expressions can be used for other type 
of systems, such as atom- molecule or molecule-molecule 
collisions assuming full inelastic or reactive loss, if the 
corresponding Cq coefficients are known. For example, 
the analytical expressions can be applied to collision of 
non-ground state molecules of NaK, NaRb, NaCs, KCs 
and RbCs. The present study provides useful informa- 
tion about collisional properties of hctcronuclear alkali 
polar molecules for which increasingly experimental in- 
terest is devoted. Future studies will consider the vibra- 
tional and rotational dependence of the Cq coefficient of 
the heteronuclear alkali molecules, the higher anisotropic 
terms in the long-range interaction, as well as the ef- 
fect of higher collision energies, when more partial waves 
dominate. 
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Appendix A: Characteristics of the heteronuclear 
alkali molecules 



We provide in Table IIVI a summary of the characteris- 
tics of the fcrmionic and bosonic isotopes studied in this 
work. Conversion factors from atomic units (a.u.) to 
S.I. units are: I a.u. of mass is equal to 1822.89 a.m.u. 
(atomic mass unit), 1 a.u. of electric dipole moment is 
equal to 2.5417 D, 1 a.u. of C§ is equal to 1 EhaQ with 
1 E h (Hartree) equal to 4.35974394x 1(T 18 J and 1 a 
(Bohr radius) equal to 0.529177x 10" 10 m. 



TABLE IV: Fermionic (F) or bosonic (B) character, isotope, 
reduced molecule-molecule mass /x (in a.u.), C@ coefficient (in 
a.u.) and permanent electric dipole moment d p (in D) for the 
different heteronuclear alkali molecules. 



F/B 


isotope 


fi (a.u.) 


C 6 (a.u.) 


d p (D) 


F 
B 


6 Li 23 Na 
7 Li 23 Na 


26436 
27349 


3880 


0.531 


F 
B 


40 K 87 Rb 
41 K 87 Rb 


115638 
116547 


16133 


0.566 


F 
B 


7 Li 40 K 
6 Li 4o K 


42820 
41907 


524000 


3.513 


F 
B 


6 Li 87 Rb 
7 Li 87 Rb 


84695 
85608 


1070000 


4.046 


F 
B 


6 Li 133 Cs 
7 Li 133 Cs 


126618 
127531 


3840000 


5.355 



Appendix B: Height of the adiabatic barrier for the 
| il/i, | — 1 component in electric field. Mixing L = 1 
and L = 3. 

In this case, we have two diabatic effective potential 
curves 



V L =i(R) 
V l=3 (R) 



and a coupling 



W(R) 



2h 2 


C 6 (2/5) d 2 


2fiR 2 


R 6 Atteq R 3 


12 h 2 


C 6 (2/5) d 2 


2fiR 2 


R 6 4ire R 3 


_ (2% 


/T26/35) d 2 




47T£ R 3 



(33) 
(34) 
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In the case of \W\ <C | Vl=3~ Vl=i|, the adiabatic effective 
potential curves are given after diagonalisation by 



E±(R) = V L=3/1 (R)± 



72 fid 4 



=3/1 v-v- 875ft 2 (47re ) 2 i? 4 
and especially the lower one 

(2/5) d 2 



E-(R) 



2h 2 



+ 



2\iR 2 R 6 Aire R 3 



9± 
r 4 



with 



72 fid 4 



875 h 2 (47re ) 2 ' 



(35) 



(36) 



(37) 



At large d, the most repulsive potential in Eq. (|36|) is 
[(2/5) d 2 ]/[4TTeoR 3 } and the most attractive is -C 4 /i? 4 
so that the height of the barrier is 



E 



n— 4,m— 3 

"l=l,|M r ,|=l 



(875/24) 3 /t 6 / d 2 V 
10000 fi 3 \4tteo) 



(38) 



Appendix C: Adiabatic potential for the \Ml\ = 
component in electric field. Mixing L — and L = 2. 

Now we have the two diabatic effective potential curves 



Vl=o(*) = 



Vb=a(i2) 



C 6 (4/7) d 2 



6 ft 2 



2^ii? 2 i? 6 47re R 3 



and the coupling between them 

2d 2 



W(R) = - 



V^ineo R 3 ' 



(39) 



(40) 



In the case of \ W\ <C | Vl=i— Vl=o|j the adiabatic effective 
potential curves are given after diagonalisation by 



E ± (R) = V L=2/0 (R)± 



A/id 4 



2 M-/- i5h 2 (47re ) 2 i?4 
and especially the lower one 



with 



\\id 4 



(41) 



(42) 



(43) 



15 h 2 (4tt£ ) 2 ' 

Appendix D: QT expression for imaginary scattering 
lengths and scattering volumes 

We provide here the analytical QT expressions for 
imaginary scattering lengths and imaginary scattering 



volumes. If we define the scattering length and the scat- 
tering volume (see Ref. J6(|) by 



a 
V 



a r — i a,i = —8(k)/k 
V r -iVi = -6{k)/k 3 , 



(44) 
(45) 



for vanishing wave- vectors k — > 0, the loss rate can be 
written as 



i3 L=0 = yAhir^/fij x A 
Pl=i,m l = ( 4: fi tt k 2 Vi / fi \ xA (46) 



for one component Ml. Similarly, the elastic rate is given 
by 



#=0 = [±hnk\a\ 2 /fi\ x A 



ILiMl = [^k 5 \V\ 2 /fi\ x A. (47) 



To get the corresponding cross sections, one has to divide 
the rates by the relative velocity v = hk/fi. Identifying 
the loss rate with the QT model, one gets the imaginary 
scattering length in the van der Waals regime 



, M 1/4 C fi 
1.92 x ' ^ 



1/4 



2 7 /4 ftl/2 



(48) 



the imaginary scattering length in the electric field 
regime 



3.74 x 



/' 



H 2 V30 J 47T£ 



(49) 



the imaginary scattering volume in the van der Waals 
regime 



Vi = 0.53 x 



/ 3 9 /y/ 4 g 6 3/4 
32 ft 3 / 2 



(50) 



and the imaginary scattering volume in the electric field 
regime 



V, = 0.54 x 



3 g /y \ d 6 

ft 6 2V2 53 ) (47re ) 



(51) 



In the case of lossy collisions, the imaginary parts a% or 
Vi contributes to the elastic part of the rates. As a con- 
sequence, they provide a minimum value for the elastic 
rates /3f, =0 = (4 h tt k a 2 / fi) x A for s-wave collisions and 
/3f-! m l — (4S.7rfc 5 V 2 /fi) x A for p-wave collisions. In 
other words, lossy collisions imply non-zero elastic cross 
sections or rate coefficients. 
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